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We generalize Ehrenfest’s equations to systems having two work terms, i.e. systems with three 
degrees of freedom. For black holes with two work terms we obtain nine equations instead of two 
to be satisfied at the critical point of a second order phase transition. We hnally generalize this 
method to a system with an arbitrary number of degrees of freedom and found there is 
equations to be satisEed at the point of a second order phase transition where N is number of work 
terms in the first law of thermodynamics. 


I. INTRODUCTION 

Black hole thermodynamics has become a very 
important topic over the past few decades [I|. As a 
thermodynamic system, black holes undergo a phase 
transition shown by Davies for the Kerr-Newman black 
hole [^. Using the path integral approach, Hawk¬ 
ing and Page investigated Schwarzschild-AdS black 
holes and found the critical tem^rature and mass at 
which a phase transition occurs (3|. This subject has 
ever since been considered by a number of authors imi. 

It has been shown that the thermal AdS space rnw 
collapse to form a black hole at certain temperatures Q ■ 
This is an example of first order phase transition [13 and 
can also occur for the case of charged and Kerr-Newman 
black holes in the AdS space [l^, [T^ . The authors of [l^ 
argued that this kind of transition may also occur in flat 
and de Sitter spaces. They also showed that for black 
holes below the critical temperature of Hawking-Page 
phase transition, or the so called ’supercooled black 
holes’, there is a line of first order phase transition in 
the phase diagram that terminates at a second order 
transition point. 

The aim of this paper is to consider the critical point 
at which heat capacity diverges. In the laboratory 
thermodynamics of the second order phase transition, 
there is a finite discontinuity in the diagram of specific 
heat. But as it has been shown by Davies 0 that an 
infinite discontinuity also exist for the critical point in 
black hole’s thermodynamics. So, one can argue that it 
is not a second order phase transition analogous to the 
ones in ordinary thermodynamics [2^. 

Discontinuity of specific heat has made it interesting 
to study the critical phenomena at the divergence point 
[il[i3,[I3- Specific heat is approximated around the 
critical point to a power of t emp erature, and then the 
critical exponent a is read of [^. In it was shown 
that for the case of charged black holes, the critical 
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exponent does not depend on the dimensions of the 
space. Specific heat is also important for the study of 
local stability of black holes in the sense that black 
hole is locally stable if its heat capacity is positive. 

Here, we use the Clapeyron’s and Ehrenfest’s scheme 
to investigate the critical point. In [l^ Clapeyron 
equation is used to study the coexistence line of first 
order phase transition. Ehrenfest’s equations must be 
satisfied at the critical point of a second order phase 
transition. The critical point is the one in which heat 
capacity is discontinuous and, as we will see, this is the 
case for a large class of black holes. 

For some years, it was thought that Ehrenfest’s 
equations are not both satisfied at the divergent point 
of heat capacity; rather numerical investigation was 
used to show that black holes undergo a glassy phase 
transition [23,123. Later, examination of the divergence 
point of other parameters that appear in Ehrenfest’s 
equations showed that there was no violation in Ehren¬ 
fest’s equations, so that the transition is of second order 
(23. In this paper, we check Ehrenfest’s equations 
in the point of phase transition in a large class of 
black holes. Also, we will obtain a generalized form of 
Ehrenfest’s equations for the case of black holes that 
have more than one work term in their associated first 
law of thermodynamics and check the relations for phase 
transition in a Kerr-Newman black hole. It is fascinating 
that all of the nine generalized Ehrenfest’s equations are 
consistently satisfied at the point of second order phase 
transition. Finally we generalized Ehrefest equations to 
the case of arbitrary number of degrees of freedom. 

This paper is organized as follows. In Section H, we 
review and check the validity of Ehrenfest’s equations for 
a large number of black holes with two degrees of free¬ 
dom. In Section HI, we generalize Ehrenfest’s equations 
for systems with three degrees of freedom and obtain nine 
equations to be satisfied at the point of second order 
phase transition. Then, we check their validity for Kerr- 
Newman black hole at the divergence point of specific 
heat in different ensembles. We generalize our approach 
to the case of N work terms in Section IV and conclude 
the paper in Section V. 
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II. SECOND ORDER PHASE TRANSITION 
FOR TWO PARAMETER BLACK HOLES 

Phase transition in thermodynamic systems is an im¬ 
portant area in physics. Ehrenfest categorized phase 
transitions by using the free energy as a function of ther¬ 
modynamic variables. In this way, the order of phase 
transitions is defined as the smallest order of free energy 
differentials that have noncontinuous behavior through 
the transition. For example, the transition of solid-liquid- 
gas in water is a first order phase transition because 
entropy and volume (which are the first differentials of 
Gibbs free energy) are noncontinuous through the transi¬ 
tion. In this Section, we investigate the thermodynamic 
phase transition in a variety of black holes that have just 
one work term in their associated first law of thermody¬ 
namics. Consider the case of charged black holes. The 
first law of thermodynamics for this case will be of the 
following form: 


dM = TdS + (1) 

in which, M is the mass of black hole, T is Hawking tem¬ 
perature, S is entropy, and $ and Q are electric potential 
and charge, respectively. The work term in this equation, 
which is analogous to —PdV in ordinary thermo¬ 
dynamics, is important for writing Ehrenfest’s equations. 
As we will show soon, Ehrenfest’s equations are necessary 
conditions for the existence of second order phase tran¬ 
sitions. To obtain these equations, we need to develop a 
suitable form of free energy which is one that ensures the 
second law of thermodynamics (i.e. increasing the total 
entropy when it decreases). 

In canonical ensemble, the suitable form of free energy is 
defined as follows: 


G = M-TS (2) 

Considering o, we have: 

dG = -SdT + (3) 


We can expand the first differentials S and $ in terms of 
independent parameters Q and T: 


dS=i^)QdT+{^)TdQ = ^dT + <^adT (4) 


where, k = From the continuity of $ in phase 

transition, it follows that (i$i = d^ 2 - So, we find the 
second Ehrenfest’s equation as follows: 




Q!2 ~ Oil 
K2 — Ki 


(7) 


For the case of rotating black holes, Q and 4) are replaced 
with angular momentum J, and angular velocity O, 
respectively. 

The order of phase transitions can be identified by 
checking the validity of Ehrenfest’s equations at the crit¬ 
ical point in which heat capacity diverges. Satisfaction 
of Ehrenfest’s equations has been shown for Reissner- 
Nordstrom and Kerr black holes in AdS space [U, [2^ . 
We check these equations for Einstein-Gauss-Bonnet 
(EGB), Einstein-Maxwell-Gauss-Bonnet {EMGB), 
Einstein-Yang-Mills-Gauss-Bonnet (EYMGB), Kerr, 
Reissner-Nordstrom {RN), Kerr-AdS, RN-AdS, and 
Horava-Lifshitz (HE) with k = —1 black holes in the 
canonical ensemble. 


Now, let us consider the grand canonical ensemble, 
in which J or Q are no longer constant and can be 
exchanged between the black hole and its environment 
while their conjugate parameters are fixed. By choos¬ 
ing the suitable form of free energy, which is G = 
M — TS — $(5, we can find Ehrenfest’s equations for 
the charged black holes in the grand canonical ensemble: 


.9$ (G$)2-(G$)i 

TQ{a2-ai) 
a2 — ai 
^ K2 - ’ 


( 8 ) 

(9) 


where, G^ is the specific heat at a constant electric po¬ 
tential, a = ^(f§)T is analogous to volume expansion, 

and K = ^(f$)T is analogous to isothermal compress¬ 
ibility. It is straightforward to prove that both Ehrenfest 
equations ([5]) and ([5]) are satisfied for EMGB, EYMGB 
and HL black holes. 


HI. GENERALIZED EHRENFEST’S 
EQUATIONS FOR THREE PARAMETER BLACK 
HOLES 


where, Gq is the specific heat in constant charge and 
«= ¥(||)t- At the point of a second order phase tran¬ 
sition, we should have continuous S and $ through the 
transition so that Si = S 2 (where the subscripts 1 and 2 
are the phases before and after transition). As a result, 
dSi = dS' 2 , and from o we will have: 


(Gq)2 - {Cq)i 

$r(a2-ai) 


(5) 


which is the first Ehrenfest’s equation. In the same way, 
we have for $: 


( 6 ) 


In this section, we will generalize the Ehrenfest equa¬ 
tions to the case of three parameter black holes that have 
two work terms in their associated first law of thermody¬ 
namics. An example of this type of black hole is Kerr- 
Newman black hole which have both electrical and an¬ 
gular terms in it’s mass formula. After a brief notation 
of thermodynamics of Kerr-Newman black hole we mod¬ 
ify Ehrenfest’s equations and we obtain nine equations 
which are supposed to be satisfied at the critical point of 
second order phase transition. The Smarr mass relation 
for the Kerr-Newman black hole states that [ 2 ^: 

M = {2S + ^{J^ + ^Q^) + ^Q^)y ( 10 ) 


d^ = —^adT + ^ndQ 
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The first law of thermodynamic is 

dM = TdS + ndJ + ^dQ 


( 11 ) 


Using the first law of thermodynamics one can determine 
temperature, T, angular velocity, U, and electrical poten¬ 
tial, $. 

The suitable form of free energy for canonical ensemble 
is: 


G = M-TS. 

So, using Equation (ED we will have: 

dG = -SdT + ndJ + 


( 12 ) 

(13) 


We see that G is a function of T, J, Q. We can write dG 
as follows: 

r)G r)G r)G 

dG = {^)j,QdT+i — )T,QdJ+{ — )TjdQ. (14) 

Comparing this equation with (fT3l) . we obtain 

^ .dG. ^ .dG. ^ .dG. 

S = -i^)j.Q, U=(—)t,q, ^ = { — 

And thereby we obtain Maxwell’s equations: 

dS 

(^)t. = -(^).q, 


dn _ 9$ 

(^)t,j - (^)t.q. 


(16) 


Ehrenfest equations can be obtained by expanding the 
first differentials of free energy {S, U and $) in terms of 
the independent parameters (T, J, and Q) and by fixing 
the differentials at the transition point. For 5, we have: 


dS = {-^)j,QdT + {-^)T,QdJ -f {-^)T,jdQ 


Cj.Q 

T 


dT + fladJ + ^a'dQ 


(17) 


where, 


l,dS. 1,5U^ 

, 1 95 19$ 

“ = TfT(A77)uJ =-;^(^)J.Q 


Since entropy is continuous, we have dSi = dS 2 , in which 
the indices show the states before and after the transition. 

(GyQ)i + ™«l(^)s + 

= (GyQ)2+raa2(^)s + T$a'(^)s (18) 

which can be written in the following form: 

__ (Cj,q)2 - (Gj,q)i $(a2-a'i) dQ 

^dT^^ nT{a2-ai) n{a2 - aiy dT ’^ ^ ^ 


The general expansion of J in terms of T, 5, and Q is: 

7 T / ^\ im / ^'^ \ in / ^^ \ m\ 

dJ = {-^)s,QdT + {-^)T,QdS + {-^)s,TdQ 

By, differentiating J with respect to T for a constant S, 
we have: 

,dJ /9J. /9J, ,dQ. . . 

Using m and (EOD, we find the first two generalized 
Ehrenfest’s equations for systems with three degrees of 
freedom: 


_,dl (Gj,q)2 - (Gj.q)i 

{qt)S,Q nT{a2-ai) 

_,d£ ^ $(a'2 - a[) 

n{a2-ai)' 


( 21 ) 

( 22 ) 


Also, from we can write the following relation in¬ 
stead of (IT^ : 

_ ^ _ (Cj,q )2 - (Gj,q)i , U(a 2 - Qi) 

^dT’^ $r(a^-ai) ^ ^{a'2-a[ydT’^^ ^ 

It is straightforward to show that: 

Comparing (1^ and (HI, we find the third and forth 
generalized Ehrenfest equations: 


{Cj,q )2 - (Gj.q)i 
{qt)S,J <^>T{a'2-a[) 
,dQ _ U(q!2 - ai) 


(25) 

(26) 


Now we write U as a function of T, J and Q, so we have: 

dVL = —fladT + ^ndJ + ^x'dQ- (27) 

Where 

1 dn , 1 9fl 


$'9Q^ 


We can repeat the above calculations for U. The result 
is four Ehrenfest equations: 


dJ a2-ai 

{gj.h.Q - 

_,dJ_ ^X2 - X'l) 

U(«:2-ft:i) 

.9g _ U(a2 - ai) 

^dT>^’^~ ^X2-x'l) 
,dQ. _ n{K2 — Ki) 

ft T — IfJTT 


(28) 

(29) 

(30) 

(31) 


$(X'2-X'i)’ 

For $ as a function of T, J and Q, we have: 

(i$ = —^a'dT -\- nxdJ + ^ndQ. (32) 
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Where 


9$ 

1 9$ 


= 

- a'l) 

(33) 

to 

1 

dJ 

- k[) 

^{X2 - Xl) 

(34) 

(dQ 

- 

a '2 - a{ 

(35) 

K '2 - 

dQ 

^(X2 - Xl) 

- k[) ■ 

(36) 


All the twelve generalized Ehrenfest’s equations are pre¬ 


sented in Table ID It should be noted that some of the 
equations are the same and that a total number of nine 
independent equations exist. 

The equations (ESI) and (1551) convert to the usual 
Ehrenfest’s equations in the limit of Reissner-Nordestrom 
black hole. Also, the equations (l2T|l and (l28ll go to the 
usual Ehrenfest’s equations in the limit of the two para¬ 
metric Kerr black hole. 

Explicit expressions of singular parameters appearing in 
the nine generalized Ehrenfest’s equations are given in 
Table IIIII A is a common factor in the denominator of all 
the singular parameters. This factor is zero at the critical 
point. Canceling out this factor is the key to the satis¬ 
faction of Ehrenfest’s equations. Now we explain how 
Equation ((5T|) is satisfied. Consider the right hand side 
of Equation (EH- Using the equations for Cyg and a 
from Table [TTIl we have: 


(Cj,q )2 - (Cj,q)i 
UT{a2 — ai) 


2S{-U^ -Q^ + 645")(4J" + (Q^ + 8^)^) ^ j; - jr 

8^8J^S^ + 2S^{Q^ + 85)2(4J2 + (Q2 + sS){Q^ + 24,S)) "" i; - 

2V2S{4J^S + Q^ + 16Q^S^ + 64S^)i dJ 

J(4J2 + Q^ + 32Q25' + 192S'2) “ 


(37) 


So the first generalized Ehrenfest’s equation is satisfied. 
In the same way, we can show that other equations in 
Table U are also satisfied. 

We can use several ensembles to study the thermody¬ 
namics of the black hole. In the grand canonical ensem¬ 
ble, the angular velocity U and the electrical potential $ 
are constants. We also have two mixed ensembles, one is 
defined by taking J and $ to be constant and the other 
is defined by demanding that Q and U be constant. We 
called them J-fixed and Q-fixed ensemble, respectively. 
To check the existence of second order phase transition 
in mixed ensembles we should find generalized Ehrenfest 
equations in these ensembles by a similar calculations to 
the previous one for canonical ensemble. But there is 
a simple method to obtain these relations by comparing 
energy equations of both ensembles. The suitable form 
of free energy for Q-fixed ensemble is: 

G = M -TS - jn (38) 

Using the first law of thermodynamics, dM = TdS + 
Q.dJ + ^dQ we have: 

dG = -SdT - JdVL + ^dQ (39) 

Comparing this relation and m we see that by replac¬ 
ing {J ——ft and U —J} or {J —fl and il —J} 
we can obtain differential form of free energy of Q-fixed 
ensemble from the canonical ensemble. We have already 
obtained the Ehrenfest equations from free energy (ta¬ 
ble m. So by the above replacements it is possible to 
obtain a new set of Ehrefest equations in Q-fixed en¬ 
semble. In a similar way by replacing {Q —>■ — and 


$ —>■ Q} or {Q —>• $ and $ —>• —Q} we can obtain an¬ 
other set of generalized Ehrefest equations in J-fixed en¬ 
semble. By imposing both of the following replacements 
({{J —>■ —fl, U —J}, {Q —>■ —$ and $ —>• Q}} or 
{{(J —U, U — —J}^ {Q —5“ ‘h and —>■ —Q}}) on 
canonical ensemble we can find other relations for the 
grand canonical ensemble. Also by replacing {Q O J 
and $ y-)- n} mixed ensembles convert to each other. So 
table U can be used to find Ehrenfest equations of all en¬ 
sembles. We can easily demonstrate that the generalized 
Ehrenfest’s equations in the Q-fixed and J-fixed mixed 
ensembles are satisfied. But as shown in , the specific 
heat for constant U and does not diverge and, thus, so 
the possibility of a second order transition in the grand 
canonical ensemble is ruled out. 


IV. HIGHER DIMENSIONAL 
GENERALIZATION OF EHRENFEST’S 
EQUATIONS 

We can generalize the method used in the previous 
section to the case of N work terms in the first law of 
thermodynamics. We assume that we have a free energy 
with N work terms as below: 

N 

dG = -SdT 4- il,dJ^ (40) 

i=l 

To obtain Ehrenfest’s equations, we should first write the 
differentials of the free energy that are S and SliS in terms 
of independent parameters (i.e., T and Jis). So, we have 













5 


TABLE I. Generalized Ehrenfest’s equations in the canonical ensemble. There are twelve equations which nine of them are 
independent. 


S fixed 

fixed 

4> fixed 

/ dJ \ (Cj,q)2-{Cj^q)i 

f 'ir, a! 2 -a!i 


(eTT,Q nT(a 2 -ai) 

(dJ) 

Q(a2-c^l) 

(dQ) n(o^2-o^i) 

/ dQ\ CC'j.q)2-(C’j,q)i 

(grpjn.Q K2-K1 

(dJ) ‘S’ixi-xi) 

[gQjn,T n(K2-«i) 

fSQ'l 0(02-01) 

'.arl-i'.Q n(x2-xi) 

/ aj \ *(k 2 -k'i) 

laol^T n(x2-xi) 

taQ\__ _o(x2-xi) 

taQ\ __ _ a' 2 -a'i 


<I>(X'2-Xi) 

{gpJ<P,J k'2-k[ 



1 /^4>^ 

^ ^ n (ej )tQ 

« - 

O’ 

h 

|<ti 

10 

11 

^ f dS \ 

1 / dQ\ 

1 / \ 


A ^ ( 80 )^>-^ 

^ ^ « 1 80 


TABLE 11. The singular parameters appearing in the Generalized Ehrenfest’s equations for the canonical ensemble. 


Cj,Q 

25(-4J^-Q‘*+645^)(4J^ + (Q^+85)^) 

.4 

a 

8x/8-f^S’*+2S3(Q2+8S)2(4j2 + (Q2+8S)(Q^+24S)) 

4 

OL 

8x/8J2S3+2S3(Q2+8S)2(4J2(q2 8S)-I-(Q2+8S)3) 

(Q2+8S)A 

K 

(4J2-(3Q2-8S)(Q2+8S))(4J2 + (q2+8S)2) 

/ 

K 

(Q^-85)(Q^+85)''+48J^(3Q^+85)+8J^(Q"'+85)(5Q'‘+48Q"'5+1925^) 

Q(Q2+8S)A 

X 

8Q(Q2+4S)(4j2_,_(q^+8S)2) 

4 

x' 

-168(0'’+4S)(4J"’ + (Q"’+8S)-') 

(Q2+8S)A 


A = 48J'‘ + (3Q^ - 8S)(Q^ + 8Sf + 8J^{3Q'^ + 320^3 + 1925'^). 


the following relations: 



We have TV + 1 equations similar to dm, differentiat¬ 
ing any of which with respect to T leads to Ehrenfest’s 
equations. We can choose any one of N work terms in 
Equations m to write the following relations: 


dJi \ 

It ) 


{fife} 


- tti -I- ^ (3ij 


dT ) 


(42) 


{fife} 


where, at and jSij are some constants. Each set of TV -|- 1 
equations of (1^ generates N equations similar to (l42ll 
each of which leads to N Ehrenfest’s equations. So far 
we have N‘^{N + 1) equations some of which are the 

same. Sentences like ( ^ ) and ( 44^ ) lead 

\OJj J J rp g 

to the same Ehrenfest’s equations under exchange of 


i O j. We have N'^{N — 1) terms like 
{i,k = 1,2,N i ^ j) and N{N — 1) terms like 
■ So, the final number of Ehrenfest’s equations 
is given by: 


N^{N + 1) - ^ {N^{N - 1) -b N{N + 1)) 

N{N + 1)^ 

2 


(43) 


It is clear that we have two equations for = 1, nine 
equations for = 2, etc. 

Also, we can exploit the same argument to determine the 
number of generalized Ehrenfest’s equations for Myers- 
Perry black holes [2^ that are the generalized form of 
Kerr-Ads black holes in higher dimensions, with rotation 
in more than one plane. In higher dimensions, the num¬ 
ber of rotating planes will be as follows: 


N < 


fd- 11 


2 
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where, d is the number of space-time dimensions and 
brackets give the integer part. For example, we have 
two rotating planes for a five dimensional space-time and 
three for a seven dimensional one. We have nine general¬ 
ized Ehrenfest’s equations for a five dimensional Myers- 
Perry (MP) black hole and 24 equations for a seven di¬ 
mensional MP black hole. 


V. CONCLUSIONS 

In this paper, we considered the thermodynamic phase 
transition in black holes strongly focusing on Ehrenfest’s 
equations. Eor a large class of black holes which have 
just one work term in their first law of thermodynamics, 
we showed that at a second order phase transition in the 
canonical and grand canonical ensemble the Ehrenfest 
equations are satisfied. 


We also generalized Ehrenfest’s equation to the case 
of systems with three degrees of freedom. We obtained 
nine Ehrenfest’s equations that are satisfied at the critical 
points. As an example, we considered the Kerr-Newman 
black hole in the canonical ensemble. We showed that the 
generalized Ehrenfest’s equations are satisfied at the crit¬ 
ical point in which heat capacity diverges. So, the infinite 
discontinuity of specific heat in black hole thermodynam¬ 
ics causes no problem for a second order phase transition 
to occur. In the case of the Kerr-Newman black hole, 
we found two mixed ensembles that have their own gen¬ 
eralized Ehrenfest’s equations. Einally, we generalized 
this approach to the case of N work terms in the first 
law of thermodynamic which can be used to check the 
order of transitions in black holes that have more than 
two work terms such as Myers-Perry black hole in higher 
dimensions. 
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